% usefile /var/folders/p1/0_bpjw1104jdbqfjqfzsxt980000gn/T//RtmpbRejU4/bayesx/bayesx.estim.input.prg

logopen using /var/folders/p1/0_bpjw1104jdbqfjqfzsxt980000gn/T//RtmpbRejU4/bayesx/bayesx.estim.input.prg.log

bayesreg b

map rcmap 
rcmap.infile using /var/folders/p1/0_bpjw1104jdbqfjqfzsxt980000gn/T//RtmpbRejU4/bayesx/rcmap.bnd
map rcmap 
rcmap.infile using /var/folders/p1/0_bpjw1104jdbqfjqfzsxt980000gn/T//RtmpbRejU4/bayesx/rcmap1.bnd
map rcmap 
rcmap.infile using /var/folders/p1/0_bpjw1104jdbqfjqfzsxt980000gn/T//RtmpbRejU4/bayesx/rcmap2.bnd
map rcmap 
rcmap.infile using /var/folders/p1/0_bpjw1104jdbqfjqfzsxt980000gn/T//RtmpbRejU4/bayesx/rcmap3.bnd
dataset d 
d.infile using /var/folders/p1/0_bpjw1104jdbqfjqfzsxt980000gn/T//RtmpbRejU4/bayesx/bayesx.estim.data.raw

b.outfile = /var/folders/p1/0_bpjw1104jdbqfjqfzsxt980000gn/T//RtmpbRejU4/bayesx/bayesx.estim

b.regress prolabor = district(spatial,map=rcmap) + unionpop*district(spatial,map=rcmap) + aapct*district(spatial,map=rcmap) + urbanpct*district(spatial,map=rcmap) + dem + laborcomm, family=binomialprobit iterations=250000 burnin=2000 step=10 setseed=998307315 predict using d 

b.getsample

logclose 
